Two pressure-induced transitions in TiOCl: Mott insulator to anisotropic metal 
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Using Car-Parrinello molecular dynamics calculations we investigate the behavior of the low- 
dimensional multiorbital Mott insulator TiOCl under pressure. We show that the system undergoes 
two consecutive phase transitions, first at P c from a Mott-insulator to a metallic phase in the ab 
plane with a strong Ti-Ti dimerization along b. At a pressure P c ' > P c the dimerization disappears 
and the system behaves as a uniform metal. This second transition has not yet been reported 
experimentally. We show that the insulator-to-metal transition at P c is driven by the widening of 
the bandwidth rather than structural changes or reduction of crystal field splittings and it shows a 
redistribution of the electronic occupation within the t2 g bands. Our computed pressure-dependent 
lattice parameters are consistent with experimental observations and the existing controversy on 
the change of crystal symmetry at high pressures is discussed. 

PACS numbers: 71.30.+h,71.15.Mb,71.15.Pd,71.70.Ch 



Understanding the Mott insulator-to-metal transition 
in strongly correlated systems is a central issue in con- 
densed matter physics due to the striking phenomena, 
like unconventional superconductivity, emerging at the 
vicinity of the transition P, 0]. Application of chem- 
ical [3, |4{ or hydrostatic pressure [5j to a Mott insu- 
lator is an ideal way to trigger the Mott transition, as 
the one-electron bandwidth (W) is directly controlled by 
pressure, leading to the change of electron correlation 
strength (U/W) with U the effective on-site Coulomb in- 
teraction. However, pressure-induced structural changes, 
as present in FeO p] or MnO Q , and orbital repopulation 
in multiorbital systems increase the degree of complex- 
ity of the transition and make it difficult to identify as a 
Mott transition. 

Recently, TiOCl, a layered Mott insulator at room 
temperature, has attracted considerable interest [1, 0, 
Hoi. [IH [l2| due to its unconventional transitions to a spin- 
Peierls state through a zero-field structural incommensu- 
rate phase. This phase has been ascribed to competing 
large magnetoelastic couplings with ferromagnetic inter- 
chain frustration [l3j |. The idea of inducing a metallic 
state or even superconductivity by application of pres- 
sure has recently been tried out by some ex per imental 
groups 14 , llj 2, • Kuntscher et al. 14 , ItJ report 
from optical measurements and X-ray diffraction analy- 
sis the observation of an insulator-to-metal transition at 
16 GPa accompanied by a structural symmetry change. 
Forthaus et al. 15], on the contrary, neither observe any 
indication of a metallic state up to pressures of 24 GPa 
from transport measurements nor do they identify from 
X-ray diffraction analysis changes on the crystal sym- 
metry at 16 GPa. Recently, Blanco-Canosa et al. [l6| 
found a dimerized insulating state at high pressure by 
X-ray diffraction and magnetization measurements com- 
bined with ab initio calculations. 

In view of the existing controversy, we present ab initio 
Car-Parrinello molecular dynamics results for TiOCl un- 



der pressure, where both lattice dynamics and electronic 
properties are treated on the same footing. Correlation 
effects are considered within the generalized gradient ap- 
proximation GGA+U functional. We find, in agreement 
with Refs.ll4l andll7l that TiOCl undergoes an insulator- 
to-metal transition at a critical pressure P = P c and the 
transition is accompanied by a lattice dimerization along 
b with a change from orthorhombic Pmmn to monoclinic 
P2i/m s ymm etry. This structural change is consistent 
with Ref.[l6|. Moreover, our results show that the system 
undergoes a second phase transition from a dimerized 
(P2i/m symmetry) to a uniform (undimcrized) metal 
(Pmmn symmetry) at a pressure P = P c ' > P c . This 
transition has not yet been reported experimentally 

The microscopic analysis of the insulator-to-metal 
transition at P = P c reveals that this transition hap- 
pens in the ab bilayers with closure of the gap at the 
Fermi level, while the conductivity along c may be hin- 
dered by the large van der Waals gap between differ- 
ent bilayers, reminiscent of the anisotropic conductivity 
in cuprates upon doping [l8j |. The transition is driven 
by the widening of the Ti 3d bands rather than a re- 
duction of the crystal field splittings and we find or- 
bital repopulation of the electronic density in the ti g 
orbitals. The concomitant dimerization along b is at- 
tributed to the pressure-enhanced hybridization. We 
show that both phase transitions at P c and P' c are of 
first order. The calculated pressure-induced changes of 
lattice parameters are in good agreement with the exper- 
imental results 3, 17] and we can explain their strong 
anisotropic properties as well as their abrupt jump be- 
havior (contraction-expansion) at the transitions. 

For our calculations, we employed the Car- 
Parrinell o [ 1 Qi] projector-augmented wave (CP-PAW) 
method [2fJ and studied the pressure effects on the 
electronic and lattice variables with the Parrinello- 
Rahman fictive Lagrangian approach 2l|. The great 



advantage of this procedure is the ability to deal with 
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FIG. 1: (Color online) Pressure dependence of lattice param- 
eters and volume. Pressure is normalized to the critical pres- 
sure for the insulator-to-metal transition and lattice parame- 
ters are normalized to their values at ambient pressure. The 
compressibility along a, b and c is strongly anisotropic. The 
existence of two phase transitions is recognizable by jumps in 
all quantities. 



crystals with anisotropic compressibilities along different 
directions as is the case for TiOCl. In order to have 
the correct insulating behavior at ambient pressure, we 
employed the GGA+U functional. It has already been 
shown that strong correlations are also important 

for the correct description of the lattice structure. In 
fact, the lattice parameters and the distances and angles 
between atoms obtained from CP-PAW with GGA+U 
at ambient pressure are more consistent with the experi- 
mental results than GGA or local density approximation 
(LDA) determined structures. We performed carefully 
converged calculations with high energy cutoffs of 45 Ry 
and 180 Ry for the wave functions and charge density 
expansion, respectively. A value of U — 1.65 eV with 
which correct spin exchange along b is reproduced (l3] | 
was used for the GGA+U calculations. The relation 
between the U value used in different codes is discussed 
in ref. 

In Fig. Q] we present the calculated change of lattice 
parameters as a function of applied external pressure. 
The lattice parameters at each pressure are obtained by 
fully relaxing the lattice structure without any symme- 
try constraints. Although a spin ordered state has to 
be imposed within the GGA+U functional, the results 
are completely consistent with the experiments [l5], [l?} • 
We observe that the largest compression happens along 
c since the interaction between neighboring bilayers is of 
van der Waals type. Along b, which is the direction of the 
Ti chains with strong antiferromagnetic (AF) interaction 
Jf, 0], a noticeable shrinkage is observed as pressure in- 
creases. The magnitude of the shrinkage depends on the 
choice of the spin configuration set in the GGA+U calcu- 
lation. The experimental results are reproduced when we 
set AF order along b while this direction becomes incom- 



pressible if the spin configuration is set to be ferromag- 
netic (FM). This is due to the fact that the Pauli principle 
prevents two electrons with the same spin to occupy the 
same bonding state with low energy. Therefore, if the 
volume along b is reduced in a FM spin configuration, 
the total energy will increase rapidly due to occupation 
of an antibonding state as two Ti ions come closer. For 
an AF configuration along b, the Pauli principle doesn't 
affect the energy and the system becomes compressible 
along b as observed [15 [ . 

In contrast, TiOCl remains almost incompressible 
along a regardless of the spin configurations we chose 
along this direction. This is due to the fact that the Ti- 
Ti interaction along a is of Tid x 2_ y 2-0 p-Tid x 2_ y 2 su- 
perexchange nature, and not direct exchange as along b. 
In this case spin order is irrelevant for compressibility. 
Throughout this work, we adopted AF order along b and 
FM along a and c 0. 

At P — P c (black vertical line in Fig. [T|) all lattice 
parameters show an abrupt jump. While b and c con- 
tract, the a direction expands, as experimental observa- 
tions seem to indicate [23(. We find a second jump of the 
lattice parameters at P = P' c > P c (blue vertical line in 
Fig. [1]) that defines a second phase transition. At this 
pressure, both a and b directions show an abrupt expan- 
sion. In order to understand the nature of these phase 
transitions, we now analyze the density of states (DOS) 
and change of symmetry. 

In Fig. [5]we show the Ti d partial DOS in the range of 
energies between -2 eV and 6 eV at pressures in the vicin- 
ity of P c and P' c . We chose the local coordinate system 
as x || b, y \\ c and z |j a, which diagonalizes the Hamilto- 
nian with d x 2_ y 2, d xz , d yz forming the ti g bands and d xy , 
d z 2 the e g bands. In the insulating phase at P/P c = 0.51 
(see Fig. [2] (a) and (b)), a gap is present in the DOS, 
separating the occupied d x 2_ y 2 state from the unoccu- 
pied d yz and d xz states. The whole character of the DOS 
and structure is similar to that at ambient pressure for 
temperatures above the spin-Peierls transition [1, [l^] ex- 
cept for a reduction of the gap and the width of the d 
bands. At P/P c = 1.03, right above the transition (see 
Fig. [2] (c) and (d)), the gap disappears and the d states 
are suddenly rearranged. The Fermi level crosses the 
d yz and d xz states while a gap within d x 2_ y 2 (orbital 
in the be plane) still remains. These results indicate an 
insulator-to-metal transition at P c . Moreover, the abrupt 
expansion along a can be explained by the occupation of 
the dyz and d xz orbitals which strongly overlap between 
neighboring Ti ions along a. Direct repulsive interac- 
tions between same-spin electrons in these orbitals based 
on the Pauli principle lead to a sudden expansion. We 
also observe a sudden contraction of the structure along 
b and a dimerized pattern, as shown in Fig.[2](d), which 
is due to the occurrence of tighter Ti-Ti dimer states be- 
tween neighboring sites along b and explains the gap in 
the d x 2_ v 2 partial DOS. The obtained insulator-to-metal 
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FIG. 2: (Color online) (a), (c), and (e) are the partial den- 
sity of states at three different values of pressure in different 
phases. In the insulating state, only d x 2_ y 2 orbital is occupied 
while all the ti a bands are occupied in the metallic regions. 
In the dimerized metallic phase, a gap still opens in d x 2_ y 2 
orbital, (b), (d), and (f) is the structure at these three differ- 
ent values of pressure. The structure in the dimerized metal 
is quite different from the other two. In the metallic phase, all 
the O atoms behind Ti atoms are invisible due to the flatness 
of the Ti-O-Ti angle. 



transition is in agreement with optical experiments 11411 
but contradicts electrical resistivity measurements |15| . 
TiOCl is a strongly anisotropic crystal with a large van 
der Waals gap along c. If the resistivity measurements 
have been done along c, the predominance of this gap is 
to be expected. We observe from our calculations that 
the metalicity is in the ab bilayers 



14] 



In Fig. [U (e) we show the DOS at P = 1.54P C right 
above the transition at P = P c '. The gap in d x 2_ y 2 
as well as the structural dimerization along b disappear 
(Fig. [2] (d)). The system shows a uniform metallic phase 
which has not yet been established experimentally. The 
abrupt increase of occupation numbers in d yz and d xz or- 
bitals at P = P' c accounts for another sudden expansion 
in a, and the disappearance of the dimer states along b 
for the abrupt expansion of b. 

In Fig. [3J we present the pressure-induced changes of 
the d-state bandwidths as well as the crystal field split- 
tings between the lowest d-band (d x 2_ y 2) and all other 
t2 g and e g bands. We observe within the three phases 
that the bandwidths of all d states widen monotonically 
with pressure (Fig. [3J (a)), while the crystal field split- 
tings show two different trends (Fig. [3J (b)). The orbital 
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FIG. 3: (Color online) (a) Bandwidth of all five d states as 
a function of applied pressure. Pressure is normalized to the 
critical pressure, bandwidth to the value at ambient pressure, 
(b) Crystal field splittings from d x 2_ y 2 to the other d orbitals 
as a function of pressure. Lower symmetry in the dimerized 
metallic phase produces two values instead of one for both 
quantities. 



excitations from d x 2_ y 2 to the e g orbitals increase mono- 
tonically with pressure, in agreement with experiment 
(for the optically allowed excitation d x 2_ y 2 to d xy ) [Til ]. 
Meanwhile, the excitations to the t^g orbitals depend 
very little on pressure. From these results we conclude 
that the insulator-to-metal transition is mainly controlled 
by an increase of bandwidths rather than by crystal field 
splitting. However, the phase transition is characterized 
by a redistribution of electronic occupation among all 
three ti g bands. Thus, it is not a typical one-band Mott 
transition where electron transfer occurs only within one 
band. Moreover, in agreement with the sudden changes 
of lattice parameters at pressures P = P c and P = P' c , 
bandwidths and crystal field splittings also show abrupt 
jumps at the phase transitions. The cooperative effect 
between variations of electronic bandwidths and changes 
of crystal field splittings by lattice distortions may be 
the mechanism responsible for the first order phase tran- 
sitions. 

The analysis of Ti-O, Ti-Cl distances as well as Ti-O-Ti 
angles as a function of pressure shows that, as expected, 
the distances between atoms in the insulating phase de- 
crease monotonically with increasing pressure. Nonethe- 
less, the distances between Ti-0 along a increase sud- 
denly at pressures P c and P c ' due to the expansion along 
this direction at the phase transitions. The strongly 
distorted octahedron [TiO^y at ambient pressure is 
modified towards a regular octahedron by a flattening 
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FIG. 4: (Color online) angle between Ti-O-Ti as a function of 
pressure normalized by the critical pressure. The inset shows 
the difference between long and short Ti-Ti bonds along b. 



of the Ti-O-Ti angle with increasing pressure (Fig. 0J). 
At P > P^, the angles become slightly larger than 180 
degree. 

A subject of controversy has been the crystal symmetry 
for pressures above P c [l4, T^, 17 1. We have identified the 
calculated structure as orthorhombic space group Pmmn 
(59) for pressures P < P c and P > P' c and monoclinic 
space group P2\jm (11) for pressures P c < P < P c '. In 
order to confirm that the restoration of the high symme- 
try Pmmn at pressures P > P' c is independent of the 
choice of unit cell and spin order, we also performed cal- 
culations with different spin configurations (AF/FM) for 
different supercells and artificially disordered structures 
at pressures P > P' c . The high symmetry Pmmn is re- 
covered in all cases. These investigations indicate that 
the lattice first undergoes a transition from orthorhom- 
bic to monoclinic structure as dimerization sets in, which 
is consistent with Ref. 0, and under further pressure in- 
crease a transition back to the orthorhombic symmetry 
occurs as the dimerization vanishes due to nonlineari- 
ties in the elastic response. The abrupt onset and dis- 
appearance of dimerization shown in the inset of Fig. [4] 
indicates that both phase transitions are of first order. 
Such complicated structural phase transitions and high 
lattice anisotropy could be the origin of the difficulties 
in determining the lattice structure from X-ray diffrac- 
tion data, leadin g to different conclusions from experi- 
ments [13, EE EEH3] ■ 

Finally, we note that the phase transitions in the CP- 
PAW calculations with the GGA+U functional occur at 
higher pressures (39 ±1 GPa) than the experimental pres- 
sures (« 16 GPa). Such deviations are attributed to the 
artificial long-range AF order along b we imposed, which 
is used to mimic the strong correlation in the GGA+U 
approach. However, as is well-known, only short-range 
AF fluctuations exist in a Mott insulator. The situation 
can be clearer if we resort to the phase diagram shown in 
Ref. [HI]- It is found that one has to apply much higher 
pressure to drive an AF insulator to a Fermi liquid than 
to drive a Mott insulator to a metal. To remedy this shift 
in the transition pressure, relaxation with a yet to be de- 



veloped constant pressure ab initio molecular dynamics 
LDA+DMFT method would be required. 

In conclusion, our ab initio molecular dynamics study 
of TiOCl under pressure shows that this system under- 
goes two phase transitions, the first one from an insula- 
tor to an anisotropic metal accompanied with a strong 
structural Ti-Ti dimerization along b, and a second one 
from the dimerized to a uniform metallic state. Our 
results indicate that the Mott insulator-to-metal tran- 
sition is mainly driven by bandwidth widening leading 
to an abrupt change of the electronic structure and or- 
bital population and inducing a structural change. We 
propose resistivity experiments contacting the ab planes 
in order to settle the question of insulator or metal as 
well as analysis at higher pressures to confirm our pre- 
diction of a second phase transition. Due to the orbital 
repopulation and the structural change, which leads to 
a dimensional crossover from quasi-one dimensional to 
quasi- two dimensional system at P > P c ', the question of 
possible superconductivity is still open. 
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